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Abstract 

An archetypical spin-glass metallic alloy, Cuo.ssMno.iy, is studied by means of an ab-initio based 
approach. First-principles calculations are employed to obtain effective chemical, strain-induced 
and magnetic exchange interactions, as well as static atomic displacements, and the interactions 
are subsequently used in thermodynamic simulations. It is shown that the calculated atomic and 
magnetic short-range order accurately reproduces the results of neutron-scattering experiments. 
In particular, it is confirmed that the alloy exhibits a tendency toward ordering and the corre- 
sponding ordered phase is revealed. The magnetic structure is represented by spin-spiral clusters 
accompanied by weaker ferromagnetic short-range correlations. The spin-glass transition temper- 
ature obtained in Monte Carlo simulations by a finite-size scaling technique, 57 K, is in reasonable 
agreement with experimental data, 78 K. 
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I. INTRODUCTION 



Spin-glass materials still keep on puzzling researchers with their unusual magnetic prop- 
erties. Since the beginning of 1970s, when Canella and Mydoshi observed a cusp in the 
temperature dependence of the magnetic susceptibility of a dilute solution of Fe atoms in 
Au, the theory of spin glasses has become a subject of non-trivial physical and mathemat- 
ical investigations.^ Three years after the initial discovery of the phenomenon, Edwards 
and Anderson^ attributed the origin of the cusp to the appearance of a spin-glass state, in 
which the moments of the magnetic atoms are frozen in a disordered, glassy structure. They 
showed that a substantially simplified model (EA model) could reproduce such a state and 
the model became a basis for analytical mean-field considerations. However, despite the in- 
disputable success of the mean-field theory to explain many of the spin-glass properties, the 
connection of the EA model to real experimentally investigated spin-glass materials, such as, 
for instance, dilute magnetic alloys of Mn and Fe in noble metals (Cu, Ag, and Au), remains 
in most cases subtle and unclear. A great deal of failure to establish a relation between the 
mean-field treatment and realistic models can be ascribed to the lack of knowledge of the 
detailed structure of real materials. 

Here, we present a theoretical investigation of the structure and spin-glass behav- 
ior of Cu-Mn alloys using a realistic parameter-free model based on first-principles cal- 
culations. There is a large amount of experimental information on the structural and 
spin-glass properties of dilute noble-metal-manganese alloys and specifically in Cu-Mn 
alloys.-'^'^i^i^'^i^iii'i^ii^'i^ii^'i^ii^'i^ This especially concerns the atomic and magnetic short- 
range orders (ASRO and MSRO) in Cuo.83Mno.17 alloys which have been intensively in- 
vestigated by various experimental techniques. In particular, it has been established that 
the ASRO in Cuo.83Mno.17 exhibits maxima at Qa = (1,1/2, 0),-^'^^^'^^'^ while the MSRO 
indicates the existence of a strong tendency toward formation of magnetic spin-density 
wave (SDW) clusters with a wave vector Qm = (I5I/2 ± 5,0), where 6 is concentration 
dependent.-*^ In addition, static atomic displacements have just recently been investigated 
by diffuse x-ray scattering in this alloy.— Nevertheless, the detailed atomic and magnetic 
structures of these alloys at the atomic level are not known. For instance, the experimentally 
observed ASRO does not allow determining unambiguously the type of atomic ordering. A 
comprehensive analysis of the most accurate up-to-date experimental data^ shows the pres- 
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ence of elements of several (1, 1/2, 0)-type superstructures: Dia, DO22, and Pt2Mo, but none 
of them can be singled out as a candidate for the low-temperature ground-state structure. 

Magnetic ordering in these alloys also exhibits quite an unusual behavior. Although the 
neutron-scattering experiments show strong short-range correlations of the SDW type, no 
magnetic ordering occurs in the CuMn alloys at low temperatures. Instead, typical traits of 
a spin glass (SG) are observed: the cusp in the linear susceptibility,— difference in the zero- 
field-cooled (ZFC) and field-cooled (FC) magnetizations,— frequency dependence of the ac 
susceptibility,— etc. All these phenomena are connected to the SG phase transition observed 
as the divergence of the non-linear susceptibility,— which is equivalent to the divergence of 
the SG-susceptibility (see Section rvT|) .- 

The mechanism of the onset of the SG behavior in metallic alloys, and in the CuMn alloy 
in particular, is still unclear. Along with fundamental problems concerning the universality 
class of a frustrated system with long-range interactions, there is still much to be understood 
in what concerns the role of magnetic clustering in a SG transition. Indeed, the observation of 
strong short-range correlations of the SDW type below and above the transition temperature 
indicates that spin freezing is accompanied by the development of magnetic clusters. There 
is no doubt that this collective behavior contributes appreciably to the magnetic dynamics 
near the transition. It is thus necessary to have a clear picture of the magnetic ordering at 
the atomic scale. 

At the same time, comprehensive theoretical investigations, in particular those based 
on first-principles calculations, of atomic and magnetic orders in Cu-rich Cu-Mn alloys are 
practically absent. The first ab initio calculations of the magnetic exchange interactions 
in dilute Cu-Mn alloys were performed by Oswald et al!^^- on the basis of Korringa-Kohn- 
Rostoker (KKR) impurity Green's function calculations. They showed that the magnetic 
exchange interactions between Mn atoms in Cu are of the antiferromagnetic type at the 
first coordination shell, but ferromagnetic at the second and third, explaining thereby the 
existence of the weak ferromagnetic SRO observed in Cu-Mn alloys. 

The type of the magnetic and atomic short range order in Cu-rich Cu-Mn alloys has 
also been investigated by Ling and co-workers,—!^ who calculated the ASRO and param- 
agnetic susceptibilities in these alloys using a so-called S'*^^^-formalism with cavity correc- 
tions within density functional theory KKR-coherent potential approximation (DFT KKR- 
CPA).— However, such a mean-filed consideration could not provide information about the 
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detailed atomic and magnetic structures at the atomic level. Besides, their calculations for 
the CussMuis alloy could reproduce neither splitting of the (1,|,0) magnetic peak nor the 
ordering tendency for the ASRO, for which they found the phase-separation (or clustering) 
behavior in direct contradiction to the experimental data. 

In this work, we investigate atomic and magnetic orderings in the Cuo.83Mno.17 alloy on 
the basis of corresponding Monte Carlo thermodynamics simulations, with chemical and 
magnetic interactions being deduced from first-principles calculations. The outline of the 
paper is as follows: Basic models, methods, and approximations, along with the description 
of our approach, are presented in Section [Tll The calculated chemical and strain-induced 
interactions are analyzed in Section lllli In Section IIVI the results for the ASRO obtained 
from Monte Carlo simulations are compared to those provided by experiments. Cooling 
down below the ordering temperature in the Monte Carlo simulations also allows us to reveal 
the underlying ordered phase of a dilute CuMn alloy. Monte Carlo magnetic simulations 
described in Section |V] provide us with a detailed picture of the magnetic correlations in the 
system. A more thorough investigation of the SG behavior is presented in Section IVll where 
we show that the system exhibits a spin-glass transition and examine the critical behavior. 
We conclude the results in Section IVIII 

II. METHODOLOGY 

Two different types of ordering are of interest in our case: atomic or chemical, describing 
relative positions of Cu and Mn atoms on the lattice; magnetic, associated with the spin 
configuration of the Mn atoms. In general, these degrees of freedom can be interconnected 
in a quite complicated way due to the strong dependence of magnetic interactions on the 
local chemical environment of magnetic atoms, and owing to the dependence of the chem- 
ical interactions on both the local and global magnetic states, as it is the case in, e.g., 
FeCr alloys.— However, such an interconnection is relatively weak in Cu-rich Cu-Mn alloys. 
Besides, the magnetic and atomic ordering effects are well separated in temperature. 

At high temperatures relevant for the atomic local ordering, where atomic diffusion is still 
possible (above 400-500 K), the alloy is in a paramagnetic state with randomly distributed 
directions of local magnetic moments on Mn atoms. The thermally induced magnetic ex- 
citations connected with the fluctuation of the direction of spin magnetic moments on Mn 



4 



atoms are several orders of magnitude faster than the atom-vacancy interchange jumps oc- 
curring during equihbration of the alloy. It is therefore possible to average out the magnetic 
degrees of freedom and obtain effective interatomic interactions to describe the alloy ther- 
modynamics at temperatures substantially higher than the spin-freezing temperature. On 
the other hand, at low temperatures - of the order of lOOK and below - where the spin-glass 
transition is observed, atomic diffusion is practically absent^ and magnetic configurations 
should be determined for the Mn atoms fixed in their positions on the lattice. Two sepa- 
rate problems are thus considered: finding a chemical alloy configuration on the underlying 
lattice and obtaining the equilibrium ensemble of magnetic configurations for the Mn atoms 
at a temperature of interest. 



A. Atomic configurational Hamiltonian and effective cluster interactions 

The chemical, or atomic, configurational Hamiltonian used in the present work is of an 
Ising type 

where spin-variables, CTj, take on values +1 or -1 if a site i is occupied by Mn or Cu atom, 
respectively. The effective cluster m-site interactions (ECIs), v}j^\, and the pair strain- 
induced (SI) interactions, Vfj, are obtained from ab initio calculations as described below. 

The ECI have been calculated by the screened generalized perturbation method (SGPM) 
(Ref. |32i) on the basis of self-consistent electronic structure calculations of random Cu- 
Mn alloys by the exact muffin-tin orbital (EMTO) method^^ within the coherent-potential 
approximation (CPA).— Since the magnetic structure stabilizes at high temperature in a 



paramagnetic state, we have used a disordered local-moment (DLM) (Ref. l35l ) magnetic 
configuration for Mn atoms, treating the Cu-Mn alloy as a three-component system, Cu- 
Mn^^-Mn-'-, where Mn^^'-'- represent spin- up and spin-down alloy components, respectively, with 
the same number of Mn^^ and Mn^ atoms distributed randomly relative to each other. 
The ECIs of a binary alloy are then determined as the corresponding averages of the 



ECI of the initial three-component alloy (see, for instance, Ref. l32l). The intersite screening 
constants for the screened Coulomb interactions contributing to the pair interactions^!^ 
have been determined by 864-atom supercell calculations of a random Cuo.83Mno.17 alloy in 



the ferromagnetic state using the locally self-consistent Green's function (LSGF) method,— 
assuming that the screening constants do not depend on the magnetic configuration of Mn 
atoms. 

Although the size mismatch of Cu and Mn atoms is rather moderate (the atomic volume 
difference of pure 7-Mn in the paramagnetic state and Cu is less than 10%^^), and local 
lattice relaxations should consequently be relatively small and give little contribution to alloy 
energies, they should be included whenever a quantitative analysis is performed, especially in 
case of relatively small chemical effective interactions. In order to take the local relaxation 
effects into account, we have calculated the strongest and most important strain-induced 
interactions at the first three coordination shells in the dilute limit of Mn in Cu using a 
supercell approach as 

T/.si = K^^) - r(^) - 2(E^'^^ - E^'^^ ) (2) 

ij ^iirel ij'jUnrel V rel unrel/' \ / 

where E^^^ and E^j^ are the total energies of supercells with a single Mn impurity and with 
a corresponding pair of Mn impurities in pure Cu. Indices "rel" and "unrel" designate 
energies for the relaxed and unrelaxed supercells (only atomic positions are allowed to relax, 
the volume and shape of the supercells are kept fixed). 

The total energies have been calculated by the projector augmented wave (PAW) 
method^"^ within the local density approximation (LDA)^>^ as it is implemented in the 
Vienna ab initio simulation package (VASP).— "^"^ The plane- wave cut-off energy was set 
to 330 eV. If a structure has been optimized, the internal structural parameters have been 
relaxed until the Hellman-Feynman forces on each atom are less than 0.001 eV/A. Two 108-, 
and 256-atom supercells have been used to check the convergence with respect to the super- 
cell size. The presented results have been obtained for the room-temperature experimental 
lattice spacing of 3.6792 A,— unless a different lattice spacing is specified. 

B. Magnetic Hamiltonian and exchange interaction parameters 

A Heisenberg-type Hamiltonian has been used in magnetic statistical thermodynamics 
simulations 

H ^ ^ JijCiCjGiGj J (3) 

ij 
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where is a unit vector in the direction of the local magnetic moment at a site i and Jij the 
magnetic exchange interaction parameters. The occupation numbers, q, take on values 1 or 
if the Mn or Cu atom occupies a site i, respectively, [occupation numbers are related to cTj 
in the Hamiltonian ([T]) as Cj = (1 — crj)/2]. The occupation numbers in the Heisenberg Monte 
Carlo simulations have been taken for an alloy configuration obtained from an atomic Ising 
Monte Carlo simulation. The magnetic exchange interaction parameters between Mn atoms 
in the alloy have been calculated using the magnetic force theorem formalism^ implemented 
within the EMTO method. The validity of the Heisenberg description has been checked by 
comparing interactions obtained in the DLM, ferromagnetic, and antiferromagnetic states, 
as well as by comparing the results of the PAW supercell calculations with the EMTO 
magnetic force theorem results, as described below. 

We have also investigated the influence of the local environment effects and ASRO on 
the magnetic exchange interactions for a Cuo,75Mno.25 alloy. Two supercells have been gen- 
erated: one representing a completely random alloy (the SRO parameters up to the tenth 
coordination shell are kept close to zero) and another one with the SRO parameters set to 
equal those obtained from an alloy Monte Carlo simulation for 300K (the SRO parameter at 
the first coordination shell, ai ~ —0.1). The electronic structure of both supercells has been 
calculated by the LSGF method. The magnetic exchange interaction parameters in these 
supercells have been calculated for different pairs of Mn atoms having different numbers 
of Cu and Mn atoms at the first coordination shell. The relative variation of the nearest- 
neighbor exchange interaction parameter Ji has been found to be less than 5% for most 
of the pairs. The relative difference between the average values of Ji for different SROs 
turns out to be even less significant. The exchange interactions at larger coordination shells 
proved to be almost insensitive to the local environment and SRO effects. We can thus 
conclude that the fixed values of the exchange interaction parameters, independent of the 
local chemical environment and global ASRO, can safely be used in magnetic Monte Carlo 
simulations. 
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III. EFFECTIVE CLUSTER INTERACTIONS IN Cuo.ssMno 17 ALLOY 
A. Chemical interactions 

The calculated, as described in the previous section, total, chemical and, strain-induced 
interactions are shown in Fig. [H The values of the bare (unscreened) chemical interactions at 
the first two coordination shells are shown to emphasize the importance of screening effects 
in the system. As one can see, the correct account of the electrostatic contribution leads to 
the change of the interaction sign at the first coordination shell. An appreciable part of the 
concentration and lattice-constant (at 17 at % Mn) dependence of the interaction at the first 
coordination shell (depicted in Figs. [2] and [3]) is due to the one-electron energy contribution. 
The values of the effective cluster interactions up to the 80th coordination shell, which have 
been considered in the statistical thermodynamic simulations, are listed in Table |Vl 

A relatively large screening contribution to the chemical interaction at the first coor- 
dination shell suggests that the value of the interaction is quite sensitive to the screening 
parameter, which can, in fact, be quite inaccurate due to the use of the atomic sphere ap- 
proximation even with the multipole moment correction taken into account. The accuracy 
of the SGPM interaction at the first coordination shell has, therefore, been checked in su- 
percell calculations. This can be easily done in the dilute limit by considering, for instance, 
a nearest-neighbor pair of Mn atoms in a large sample of Cu. The total interaction energy 
related to the effective pair interaction in this case is 

yMn-Mn _ ^ ^^Mn-Mn _ 2^^" _|_ (4) 

where E^""^^ is the total energy of a supercell with two Mn atoms (in this particular 
case separated by one coordination shell), the total energy of the same supercell but 
containing only a single Mn atom, and E'^'^ the total energy of pure Cu (normalized on the 
same supercell). The prefactor 1/4 is due to the definition of the Ising Hamiltonian ([1]). 

In a magnetic system, Eq. (j4]) leads to two different types of interactions, depending on 
the relative orientation of the magnetic moments of the Mn atoms: ferromagnetic (FM), 
Kot'^^~^'^^5 and antiferromagnetic (AFM), The interactions are calculated from 

the total energy of the ferromagnetic, £;^^'iT-MnT^ qj, antiferromagnetic, £;MnT-Mni^ state, 
respectively. However, the interaction y*^^)"*^^^ is calculated in the DLM magnetic state, 
and in order to compare this interaction with V^^'^"^'^, proper averaging of the latter over all 



possible spin configurations of tlie pair must be performed. The interaction in the DLM state, 
relevant to our consideration at temperatures above the spin-glass transition, is therefore 
given by the following average:—*^ 

T/Mn-Mn _ / T^MnT-Mni , T^MnT-MntA fc\ 
•^tot-DLM — 2 V °* J ' 

At the same time, the difference of these two interactions is the magnetic exchange interac- 
tion parameter of the Heisenberg magnetic Hamiltonian ([3]) at the corresponding coordina- 
tion shell 

T _ T/MnT-MnJ, T^MnT~Mnt ff.\ 
•Jij — ^tot ^tot ) \^) 

where i and j are the position indices of the Mn atoms in the lattice. 
Besides, one should bear in mind that the interaction energy 

^Mn-Mn ^^^1^^^^ ^11 the 

multi-site contributions present in the system and they should be subtracted before a com- 
parison with the chemical part of the effective pair interaction is carried out 

= Ct'^-^" - \^'^''''\ (7) 

where fi™"'*' is the multi-site contribution corresponding to the right-hand side of Eq. (jll). 
It can be obtained using the GPM multi-site interactions, calculated in the dilute limit; 
the multi-site interactions are usually sufficiently accurately reproduced by the GPM since 
they do not contain an electrostatic contribution. Although the multi-site interactions are 
relatively weak (some of the strongest three- and four-site interactions for the Cuo.83Mno.17 
alloy are listed in Tables IVII and IVlip they contribute appreciably to l]™*^^*' as one can see 
in Table [H 

TABLE I: Interaction energies for a pair of Mn atoms in the dilute limit obtained from the GPM 
(y(2)-GPM|-^j^^ _^ Q-j-j from supercell calculations, Vsc^. All energies are given in mRy. The 
parameters y^^^^"^ are calculated according to Eq. 1^; O™*^^*' is the multi-site contribution. 



y(2)-GPM 


^(2) 
^ sc 


yMnT-Mn-L 


yMriT-MnT 


Ji 


j^multi 


1.085 


1.105 


-0.728 


3.040 


-3.77 


0.220 



The total energies of 256-atom supercells, ^MnT-Mnt^ ^MnT-Mni^ ^Mn^ ^^^^ ^Cu ^^^^^ ^^^^ 
calculated with the PAW method (calculation details are given in section Hi A|) for the lattice 



spacing of 3.6792 A, as has already been mentioned. These results, together with the SGPM 



interaction at the first coordination shell obtained for the same lattice spacing in the dilute 
limit of Mn, are presented in Table [B As one can see, the agreement between the SGPM 
and PAW total-energy calculations for the nearest-neighbor effective pair interaction is very 
good. Below, we come back once more to the accuracy of the SGPM interactions used in 
the present statistical thermodynamic simulations. 



B. Strain-induced interactions 



To take the effect of local lattice relaxations into account, we have calculated strain- 
induced interactions at the first three coordination shells using Eq. ([2]), as described in 
Section III A[ The interactions have been calculated in the dilute limit using the same 108- 
and 256-atom supercells as in the previous subsection. The results are presented in Table HI] 
together with the strain-induced interactions obtained by the Krivoglaz-Khachaturyan (KK) 
method^"^"^ using the experimental value for the lattice concentration expansion, u = 
~ 0.0883, for the Cun.s.'^Mnn.i? random alloy. One can see that the KK method, as it is 



implemented in Ref. |50|, quite significantly overestimates the interactions especially at the 
third and fourth coordination shells. Therefore, in the calculations of the total effective pair 
interactions we have restricted ourselves only by inclusion the strain-induced interactions 
obtained in the first-principles calculations. Thus, we assume that the contribution from 
more distant strain- induced interactions to the configuration energetics is insignificant. Let 
us note that the long-range tail of the strain-induced interactions given by the KK model 
should be overestimated because it is valid only in the dilute limit. In random alloys, on the 
other hand, the long-range tail of strain-induced interactions should actually be damped by 
fluctuations of local lattice displacements due to the corresponding fluctuations of the local 
chemical environment of individual atoms. 

Unfortunately, there exist no reliable way of getting distant strain-induced interactions 
in a random alloy from ab initio calculations. Nevertheless, strain-induced interactions 
are actually directly connected to the corresponding static lattice displacements within the 
Krivoglaz-Khachaturyan formalism. Thus, static lattice displacements in random alloys 
carry important information about strain-induced interactions. To investigate static dis- 
placements in the alloy, we have used the PAW method to calculate the relaxed atomic 
positions in a 108-atom supercell containing 90 Cu and 18 Mn atoms (corresponding to 
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TABLE II: Strain-induced interactions (in mRy) obtained from the Krivoglaz-Khachaturyan (KK) 
formalism and directly from supercell calculations. " CS" stands for the coordination shell number, 



also given by a. The parameters Va'^^ are obtained according to the procedure described in 



Ref. 



50; 3-^6 calculated using Eq. ([2]) 
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—0.100 
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-0.166 
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0.034 








7 


0.034 
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-0.074 








9 


0.091 








10 


-0.063 









16.7 at %). The chemical and magnetic configurations have been set up in such a way as 
to minimize the short-range order parameters of the ternary alloy Mn1^-Mn''--Cu. A static 
displacement for a given coordination shell has been obtained as an average over a symmetry 
group corresponding to the coordination shell. 

The calculated atomic displacements {ximn) are presented in Table llllt where they are 
compared to the recent experimental data by Schonfeld et al.^ for Cuo.83Mno.17. We also 
show the results for {xf^^'^), {xf^^^), and {x^^^) obtained from the supercell calculations 
of a single Mn (256-atom supercell) and a pair of Mn atoms (108-atom supercell), respec- 
tively. One can see that although there is a reasonable qualitative agreement between all 
the calculated results and experimental data, the quantitative differences in some cases are 
quite large. In particular, the theory and experiment predict different signs for the atomic 
displacements of Mn atoms in some cases. One can also note that Cu-Mn atomic displace- 
ments calculated for the alloy supercell and for the Mn impurity are quite different. Besides, 
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in case of impurity they decay much slower than in case of alloy. 



TABLE III: Static displacements (multiplied by 10^ for clarity) in units of the lattice constant ob- 
tained from our supercell calculations and taken from the experimental study where Georgopoulos- 
Cohen (GC) and 3A separation methods were used (see Ref. 



21 



for details). The results of the 



calculations for Mn impurity and a pair of Mn atoms at the first three coordination shells for the 
values of {xf^^^), and {x^^'^), respectively, are given in parentheses. 
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-0.12 


0.137 
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0.14 
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-0.57 


4.02 


0.83 (2.61) 


-0.83 


-0.70 


-11.82 


-6.07 


2.82 


2.10 



We would like to note that such a comparison should be made with caution because the 
magnetic states in all these cases are different: the experiment is done in the paramagnetic 
state, the impurity calculations effectively produce results for the ferromagnetic state, and 
the alloy supercell calculations are performed in a kind of a "quasi-random" magnetic state. 
The latter may represent the DLM (paramagnetic) state only if the proper averaging over 
all possible magnetic configurations is performed. In a single calculation, any pair of Mn 
atoms is either ferromagnetically or antiferromagnetically aligned rather than being in an 
average DLM state. In our view, a quite broad spectrum of the results for static atomic 
displacements reflects the complexity of the problem that requires further investigation. 



IV. ATOMIC CONFIGURATION OF CuossMno 17 ALLOY 

The atomic configuration of the Cuo.83Mno.17 alloy has been obtained in Monte Carlo 
calculations using the Metropolis algorithm in the canonical ensemble based on the Hamil- 
tonian ([1]) at 500 K, which is the annealing temperature in the diffuse scattering experiments 
Simulation boxes of different sizes (up to 32 x 32 x 32 fee elementary 



described in Ref. 
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cubic cells) have been taken to exclude finite-size effects. A run has consisted of 15000 
sweeps, with 10000 sweeps used for collecting statistics. The short-range order parameters 



•their Fourier transform is presented 
8|. Clearly, the ASRO is reproduced 



have been determined up to the 80th coordination shel 
in Fig. m along with the experimental data from Ref. 
very well up to the width of the peaks. A slight difference in the shape can be either due 
to a small error originating from the cut-off of the interaction range, and multi-site terms 
for more than 4 sites, or due to an incomplete subtraction of the magnetic contribution 
from neutron-scattering intensity data. The latter claim is supported by subtracting the 
calculated normalized atomic intensity pattern from the experimental one and noting that 
their difference matches the magnetic intensity pattern quite well. 

The ASRO in the Cuo.83Mno.17 alloy is characterized by Qa = (1, 1/2,0) type of atomic 
ordering. The local maximum of the Fourier transform of the SRO parameters at (0, 0, 0) 
is also present, yet it is very weak, implying a strong tendency toward ordering rather 
than segregation. There has been a controversy as to what kind of an ordered phase this 
ASRO represents.— "^"^"^ None of the previously proposed ordered structures (e.g., DO22, 
A2B2, Dla, etc.) satisfies all the necessary conditions, such as the absence of additional 
reflections at points inequivalent to Qa (such as for DO22) or the abundance of Cu in the 
ordered precipitates (this disfavors A2B2). To find the sought-for ordered structure, we 
have performed Monte Carlo simulations down to low temperatures. Below the transition 
temperature (Tc ~ 2507^^), an ordered phase with 25 at% Mn has precipitated. The structure 
of the phase (first proposed in Ref. |54| as an ordered phase for the NiaMo alloy) is shown in 
Fig. [5l It consists of stripes similar to the ^2-82 structure, but with half of the Mn stripes 
being replaced by Cu. 

It is worth noting, however, that the ordering energy of the obtained CuaMn phase 
(referred to as DOqq hereafter) is very close to that of the Dla (Ni4Mo) phase. It turns out 
that three-site chemical interactions play a crucial role in stabilizing the DOqq structure. 
Unfortunately, there is no simple way to check the stability of the CusMn and Cu4Mn phases 
with respect to each other, since this would require calculating the free energy for the alloy 
at different concentrations, which is a formidable task due to a strong effect of magnetism 
on the alloy thermodynamics. 

The ordering energies of the obtained DOqq phase, as well as the energies of some other 
ordered structures for the same composition, have been calculated and compared to direct 
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TABLE IV: Ordering energies of different CusMn structures: E'tot are calculated as the difference 
between the total energy of an ordered phase and a random alloy, with the Mn atoms being in the 
DLM state; i^sGPM are obtained from the effective interactions y^™-) (m = 2,3,4). 

Structure Esqpm (mRy) Etot (mRy) Esgpm - Eggg°i (mRy) Etot - Egf^° (mRy) 

Li2 1.68 2.6 2.51 2.89 

DO22 -0.83 0.26 0.54 0.54 

DOeo -1.68 -0.29 0.0 0.0 



total-energy calculations performed with the EMTO method. The results are presented in 
Table IIVI It is obvious that among the other CusMn phases, the DOqq phase is the most 
stable one at the given composition, and the energy difference with respect to the DO22 
structure ensures that minor corrections (static displacements, phonon contribution, etc.) 
cannot change the observed behavior. 

V. MAGNETIC STRUCTURE 

As a rule, accurate numerical studies of SG systems by Monte Carlo methods are com- 
putationally very demanding. However, a general picture of magnetic correlations turns out 
to be relatively insensitive to the quality of sampling. We have employed the Metropolis 
algorithm for the Hiesenberg Monte Carlo simulations to obtain pair spin-spin correlations 
in the temperature range between 300 and 20K. The atomic configurations have been ob- 
tained in the Ising atomic Monte Carlo simulations at 500 K and kept fixed during magnetic 
simulations because the atomic diffusion is practically absent at ambient and lower temper- 
atures. 

Despite frustration that impedes equilibration in simulations even above the SG tran- 
sition temperature, the spin-spin correlations tend to saturate to their equilibrium values 
much faster than quantities characterizing the SG phase. To put this another way, the au- 
tocorrelation time associated with the spin-spin correlations changes smoothly and remains 
relatively small down to the lowest temperatures owing to the absence of the second-order 
transition. We can thus make do with a simplified Monte Carlo technique to achieve a 
reasonable accuracy with a moderate number of time steps in this case. 
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To characterize a MSRO, a spin-spin correlation function C{R) = j^J^ii^i^i+R) has 
been calculated. The Fourier transform of the spin-spin correlation function obtained from 
the simulations at 20K is presented in Fig. [6] along with the experimental data on the 
corresponding polarized neutron-scattering intensities. The maxima at Qg = (1, 1/2 ± 5, 0) 
(5 ^ 0.21 for CMn = 0.172) correspond to incommensurate SDW clusters. It has already been 
contemplated^ that the local magnetic order represents a non-collinear single-Q spin-spiral 
structure, with the directions of propagation of SDWs varying in different domains. This 
type of the local magnetic order is observed in simulations at low temperatures (20K). Since 
the correlations at the first several coordination shells are strong, these clusters, rather than 
individual spins, determine the spin dynamics close to the freezing temperature T^g- Also, 
non-collinear MSRO clusters may be responsible for the chirality ordering that accompanies 
the SG transition.— 

It is worth noting that we neglect anisotropy for the pure Cuo.ssMno.ir alloy because the 
anisotropy energy for this alloy is too small (of the order of 10~^ meV per Mn atom^^) to 
have any considerable effect on the MSRO. A situation, however, can be different in CuMn 
alloys containing heavy impurities. For instance, adding even a small amount of Pt to the 
alloy results in the anisotropy energy of up to 0.1 meV per Mn atom.— 

VI. SPIN-GLASS BEHAVIOR 

In simulating a SG system, one has to take into account the fact that most of the quantities 
specific to a SG state, such as SG susceptibility, overlap parameters, etc., are quite sensitive 
to a realization of disorder- and a configurational averaging must therefore be carried out. 
The latter requires a large number of independent runs. Moreover, to obtain unbiased 
values of ensemble averages for each configuration, it is important to have samples in the 
equilibrium state. However, due to magnetic frustrations in the alloy it is hard to achieve 
the equilibration, and improved methods must be employed to get results on a reasonable 
time scale. 

Slow equilibration imposes a restriction on a system size in simulations and therefore 
finite-size effects become an important issue. To overcome the size and boundary depen- 
dences, we resort to a technique that has become a standard tool in investigations of critical 
phenomena, namely, the finite-size scaling procedure.—"^ The general idea is to obtain a 
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dimensionless correlation length, ^/L, for different system linear sizes, L (defined in terms 
of the number of elementary cubic cells of the fee lattice), and use the scaling law 



(8) 



to determine a transition temperature Tc by taking advantage of a simple corollary of the 
fact that the dimensionless correlation length becomes independent of L at T^. To calculate 
the correlation length, two independent replicas with identical atomic distributions have 
been simulated in parallel and the correlation length has been evaluated according to the 
following well-known formula:— 



1 



Xsg(O) 



- 1, 



where 



2 sin (A;^i^/2) \j xsg{K 
27r/L(l,0,0) and the SG correlation function xsG(k) is defined as 



(9) 



XsG(k) = N 



(10) 



(11) 



with indices (1) and (2) designating quantities related to the two replicas; summations 
are performed over all magnetic atoms in the alloy, (...) and [. ..]av stand for thermal 
and configurational averaging, respectively, and N = 4L^. The FSS scaling of the SG 
susceptibility xsg = Xsg(O) is given by the relation 



(12) 



Monte Carlo simulations have been performed using the heat-bath algorithm.^ To make 



calculations practically affordable, only the exchange interaction parameters with |R. 



< 



ro = 4 elementary cubic cells have been taken into account. To ensure thermalization, 
logarithmic binning has been applied to the spin-glass susceptibility Xsg(O) and a consistency 
check of the obtained configurational averages has been made in a way similar to that in 



Ref. 



63l . namely, the following two quantities have been calculated: 



Xo ^ 



(13) 



16 



and 




av 



(14) 



where t designates the time in terms of Monte Carlo steps, Iq is the equihbration time, and 
T the size of the time window during which measurements are carried out. A calculation is 
considered converged with respect to the configurational averaging if a condition Xo = Xd 
is satisfied with a high accuracy (to within 1%). To accelerate equilibration, we have used 
the overrelaxation method,— which consists in performing microcanonical steps; that is, all 
spins are flipped sequentially according to the following energy-conserving transformation: 



L sequential microcanonical steps for each heat-bath sweep have been found to be optimal 
in terms of the actual calculation time. 

The calculated size dependence of the SG susceptibility (Fig. [7]) suggests the divergent 
behavior below lOOK. More detailed information can be extracted from the size scaling 
of the dependence of the correlation length on temperature (Fig. [8]). The results show 
unequivocally that the system undergoes a transition at a finite temperature. The absence 
of a long-range order has been checked by observing the development of the correlation 
length corresponding to the SDW ordering. Although the SG transition is clearly observed, 
one cannot rule out a crossover to a marginal behavior at larger supercell sizes. The data 
exhibit a universal scaling with Tc = 57 ± 5K and u = 0.95 ±0.1 (Fig. [H]). Using these 
parameters, the SG susceptibility can also be fitted providing a parameter rj = 0.25 ± 0.1 
(Fig. [Toll . From the relation 7 = z/(2 — 77), one can also find that 7 ~ 1.7 ± 0.3. In all plots, 
the deviation of the data for L = 16 is due to the systematic underestimation of the SG 
susceptibility at low temperatures. 

An alternative (extended) scaling relation based on the high-temperature series expansion 
was proposed by Campbell et al— Unlike the conventional scaling [Eq. ([8])], the extended 
scaling applies not only to a close vicinity of the critical region but must also hold for all 
temperatures above the transition. The extended relation, that can be written as^ 




(15) 



(16) 
(17) 
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has been applied to the calculated data. The results shown in Figs. [TT] and [12] demonstrate a 
rather good scaling, although a larger temperature interval and higher accuracy are required 
to make a more thorough comparison of the two types of scalings. 

The values of the critical exponents differ from those obtained in experimental 
measurements^^ from which we have u^xp = 1.3 ± 0.15 and 7exp = 2.3 ± 0.2. One of the 
possible causes of this disagreement is a rather high sensitivity of the critical exponents to 
the scaling corrections which might be non-negligible in our calculations because the range 
of the interactions is comparable to the size of the simulation box. Another effect, which 
can strongly influence the critical behavior, is anisotropy. It is well known that even a small 
amount of anisotropy present in a real system can change the critical behavior close to the 
transition temperature, giving rise to the values of critical exponents different from those for 
an isotropic system. Unfortunately, it is practically impossible to detect any crossover from 
an isotropic to anisotropic critical behavior in small systems such as the one we have used 
here and with such a small anisotropy energy typical for pure CuMn alloys. Simulations for 
systems with an appreciable amount of anisotropy (such as, e.g., CuMna;Pty) can therefore 
be an interesting problem for future work. 

The obtained transition temperature, 57K, turns out be slightly lower than the experi- 
mental value ( 78K). One of the possible explanations for this discrepancy is the cut-off of 
the interaction range that we had to introduce. On the other hand, Monte Carlo dynam- 
ics (only local updates without overrelaxation) has revealed a strong critical slowing down 
at about a temperature equal to 80K, and this value only slightly depends on the cut-off. 
This could indicate that while dynamics is driven primarily by short-range interactions, 
the SG transition itself is very sensitive to long-range interactions which may play a rather 
important role in the onset of the SG phase in the CuMn alloys. 

In spite of a possible strong influence of long-range interactions on the SG dynamics, it 
seems unlikely that the dominating contribution to the onset of the transition comes from 
effectively infinite-range interactions because a natural cut-off can be considered originating 
from various kinds of imperfections (vacancies, impurities, static displacements, etc.) present 
in any real material but discarded in the current study. For instance, we have neglected static 
atomic displacements in the magnetic part of the problem. Although their effect should 
be very small on the effective interactions at nearest-neighbor interactions, more distant 
interaction will be definitely damped because of the corresponding exponential damping 
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of the Cu s-like states. This means that such a system will doubtfully cross over to a 
mean-field-like transition in a realistic model. 



VII. CONCLUDING REMARKS 

The basic result of this work is the structure of the CuMn alloy obtained from first- 
principles calculations. It has been shown that both the atomic and magnetic thermody- 
namics simulations accurately reproduce the corresponding experimental data. Interestingly, 
it has already been pointed out in some studies^^, that the most probable stoichiometry of 
the ordered phase of a Cu-rich CuMn alloy is Cu-25 at% Mn; however, the scattering pattern 
of known structures with this composition {DO22, L12) is not compatible with the observed 
one, and the ^42-82 structure, that lacks this drawback, has been considered as a best de- 
scription of the ASRO. In this respect, the structure DOqq, on one hand, contains 25 at% 
Mn and, on the other hand, resembles the ^2-82 structure closely and produces a similar 
scattering intensity picture with the dominant peaks at Qa = (1, 1/2,0). This reconciles a 
seeming discrepancy between the observed stoichiometry and the structure of the ASRO. 

In spite of some limitations of the approach used, a large class of metallic spin-glass mate- 
rials, e.g., AgMn, AuMn, PtMn, AuFe, etc., can be investigated in a similar fashion. Other 
contributions to the magnetic interactions might be required to be taken into consideration, 
though. Although in the CuMn alloy relativistic effects can be neglected, in systems with 
heavy elements (Au, Pt) they are important and can lead to strongly anisotropic interac- 
tions. Most of the anisotropic contributions can be treated within the GPM formalism in 



the framework of a fully relativistic code (see, e.g., Ref . l66l) . extending thus the approach to 
highly anisotropic heavy-element metallic alloys. 
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APPENDIX 

Effective cluster interactions for two-site, three-site and four-site clusters are presented 
in Tables IVTlVIII In all cases, an m-site cluster is represented by the first atom at the origin, 
(0,0,0), and (m — 1) vectors corresponding to the other atoms in the cluster as given in 
the tables. The interactions are given in millirydberg (mRy). Note that the accuracy of the 
EMTO calculations has been set to 10"'^ mRy. 
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FIG. 1: (Color online) Chemical and total (chemical and strain- induced) interactions for the first 
20 coordination shells, y^^^ and y(2).t'are ^j^g chemical interactions with and without (bare) 
the screening contribution. The total pair interaction in the disordered phase is given by the sum 

y(2) _^ ysi^ 



0.5 




^0 5 10 15 20 

Coordination shell 



23 



TABLE V: Pair effective and magnetic excfiange interactions for the first 80 coordination spheres. 



N 


(Rx , Ry , 


Rz) 


T/(2) 
(«n,fl) 


t(xc) 
(«n,fl) 


N 


{Rx 




Rz) 




j{xc) 


1 


(0.5, 


0.5, 


0.0) 


0.45035 


-2.51559 


41 


(2.5, 


2.5, 


2.0) 


0.00003 


-0.00029 


2 


(1.0, 


0.0, 


0.0) 


-0.31257 


0.83820 


42 


(3.5, 


2.0, 


0.5) 


0.00358 


-0.00001 


3 


(1.0, 


0.5, 


0.5) 


-0.20378 


0.18558 


43 


(4.0, 


0.5, 


0.5) 


0.00030 


0.00003 


4 


(1.0, 


1.0, 


0.0) 


0.27218 


-0.29230 


44 


(4.0, 


1.0, 


0.0) 


0.00038 


-0.00118 


5 


(1.5, 


0.5, 


0.0) 


0.00205 


-0.05987 


45 


(3.0, 


2.0, 


2.0) 


0.00143 


-0.00052 


6 


(1.0, 


1.0, 


1.0) 


0.04462 


-0.00211 


46 


(3.0, 


2.5, 


1.5) 


-0.00003 


0.00016 


7 


(1.5, 


1.0, 


0.5) 


-0.03562 


0.06212 


47 


(3.0, 


3.0, 


0.0) 


0.00358 


-0.01198 


8 


(2.0, 


0.0, 


0.0) 


-0.03068 


0.07661 


48 


(4.0, 


1.0, 


1.0) 


-0.00122 


-0.00056 


9 


(2.0, 


0.5, 


0.5) 


0.01150 


0.00673 


49 


(3.5, 


2.0, 


1.5) 


-0.00050 


-0.00105 


10 


(1.5, 


1.5, 


0.0) 


0.08350 


-0.06621 


50 


(4.0, 


1.5, 


0.5) 


-0.00057 


-0.00273 


11 


(2.0, 


1.0, 


0.0) 


0.01215 


-0.07616 


51 


(3.5, 


2.5, 


0.0) 


-0.00095 


-0.01165 


12 


(1.5, 


1.5, 


1.0) 


0.00675 


-0.01083 


52 


(3.0, 


3.0, 


1.0) 


0.00082 


-0.00236 


13 


(2.0, 


1.0, 


1.0) 


-0.00328 


0.01124 


53 


(3.5, 


2.5, 


1.0) 


-0.00052 


-0.00147 


14 


(2.5, 


0.5, 


0.0) 


-0.00948 


0.01708 


54 


(4.0, 


2.0, 


0.0) 
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0.00909 
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0.04217 


55 


(4.5, 


0.5, 


0.0) 


0.00025 
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-0.00019 
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(2.5, 


2.5, 


1.0) 


-0.00178 


-0.00397 


74 


(3.5, 


3.5, 


1.0) 


0.00003 


-0.00118 


35 


(3.0, 


1.5, 


1.5) 


-0.00155 


0.00064 


75 


(4.0, 


3.0, 


1.0) 


-0.00032 


-0.00073 


36 


(3.0, 


2.0, 


1.0) 


-0.00100 


-0.00282 


76 


(5.0, 


1.0, 


0.0) 


-0.00010 


-0.00014 


37 


(3.5, 


1.5, 


0.0) 


-0.00553 


0.01155 


77 


(4.5, 


2.5, 


0.0) 


0.00010 


0.00669 


38 


(3.5, 


1.5, 


1.0) 


-0.00343 


0.00050 


78 


(4.5, 


2.0, 


1.5) 


-0.00217 


0.00167 


39 


(3.0, 


2.5, 


0.5) 


0.00185 


0.01314 


79 


(5.0, 


1.0, 


1.0) 


0.00030 


0.00041 


40 


(4.0, 


0.0, 


0.0) 


0.00000 


-0.00067 


80 


(3.0, 


3.0, 


3.0) 


-0.00005 


0.00010 
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TABLE VI: Selected three-site effective interactions. The corresponding three-site clusters are 
defined by three vectors: Rq = (0, 0, 0), i?2- 



N 




R\ 

-'- ' i 






Ri 






N 




R\ 






^2 






1 


( 0-5) 


0.0, 


0.5) 


( 0-5> 


0.5, 


0.0) 


-0.01127 


15 


[ 0-5, 


0.5, 


0.0) 


( 1 


1 5 


0.5) 


00667 


2 


( 0.5, 


0.5, 


0.0) 


( 1 


0.0, 


0.0) 


-0 02395 


16 


[ 0.5, 


0.5, 


0.0) 


( 1 5 


0.5, 


1 0) 


00561 


3 


( 0.5, 


0.5, 


0.0) 


( 1 


0.5, 


0.5) 


0.04451 


17 


[ 0.5, 


0.5, 


0.0) 


(-0.5, 


1 

-L.\J, 


1 5) 


00213 


4 


( 0-5, 


0.5, 


0.0) 


( 0-5> 


0.5, 


1.0) 


-0 00140 


18 


' 1 


0.0, 


0.0) 


f 1 5 


1 


0.5) 


00135 

V/ • \J\J 


5 


( 0-5, 


0.5, 


0.0) 


(-0.5, 


1 


0.5) 


00470 


19 


[ 0.5, 


0.5, 


0.0) 


( 2 


0.0, 


0.0) 


00235 


6 


f 1 


0.0, 


0.0) 


( 0-5> 


0.5, 


1.0) 


00125 


20 


[ 0.5, 


0.5, 


0.0) 


( 1 5 


1 5 


0.0) 


-0 01867 


7 


( 0.5, 


0.5, 


1 0) 


( 1 


-0.5, 


0.5) 


-0 00014 


21 


' 1 


1 


0.0) 


( 2 


2 


0.0) 


-0 00609 


8 


( 0.5, 


0.5, 


0.0) 


( 1 


1 


0.0) 


-0.04427 


22 


[ 0.5, 


0.5, 


0.0) 


( 2 5 


2 5 


0.0) 


—0 00231 


9 


( 0.5, 


0.5, 


0.0) 


( 1-0, 


0.0, 


1.0) 


-0.00305 


23 


[ 0.5, 


0.5, 


0.0) 


(2.0, 


2.0, 


0.0) 


-0.00724 


10 


( 1-0, 


0.0, 


0.0) 


( 1-0, 


1.0, 


0.0) 


-0.00303 


24 


[ 1-0, 


1.0, 


0.0) 


(2.5, 


2.5, 


0.0) 


-0.00230 


11 


( 0.5, 


0.5, 


0.0) 


( 1-5, 


0.0, 


0.5) 


-0.00210 


25 


[ 0.5, 


0.5, 


0.0) 


( 3.0, 


3.0, 


0.0) 


-0.00055 


12 


( 0.5, 


0.5, 


0.0) 


( 1-5, 


-0.5, 


0.0) 


-0.00181 


26 


[ 1.0, 


1.0, 


0.0) 


( 3.0, 


3.0, 


0.0) 


-0.00077 


13 


( 0.5, 


0.5, 


0.0) 


( 0.5, 


0.0, 


1.5) 


-0.00324 


27 


[ 1-5, 


1.5, 


0.0) 


( 3.0, 


3.0, 


0.0) 


-0.00088 


14 


( 0.5, 


0.5, 


0.0) 


( 1-0, 


1.0, 


1.0) 


-0.00687 


28 


[ 1-5, 


1.5, 


0.0) 


( 3.0, 


3.0, 


0.0) 


0.00000 
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TABLE VII: Selected four-site effective interactions. The corresponding four-site clusters are de- 
fined by four vectors: Rq = (0,0,0), Ri, R2, R3. 



N 




Ri 




R2 






Ri 


(A) 

V ' 

(ii'lj,ii'l,ii'2,ii';s) 


1 


( 0.5, 


0.5, 


0.0) 


( 0.5, 0.0, 


0.5) 


( 0.0, 


0.5, 0.5) 


0.00418 


2 


( 0.5, 


-0.5, 


0.0) 


( 0.5, 0.5, 


0.0) 


( 0.5, 


0.0, 0.5) 


0.00459 


3 


( 0.5, 


0.5, 


0.0) 


( 1.0, 0.0, 


0.0) 


( 0.5, 


-0.5, 0.0) 


0.00164 


4 


( 0.5, 


0.0, 


0.5) 


( 1.0, 0.5, 


0.5) 


( 0.5, 


0.5, 0.0) 


0.00003 


5 


( 1-0, 


0.0, 


0.0) 


( 1.0, 0.5, 


0.5) 


( 0.5, 


0.5, 0.0) 


-0.00113 


6 


( 0.5, 


0.5, 


0.0) 


( 0.5, 1.0, 


0.5) 


( 1-0, 


0.5, 0.5) 


0.00064 


7 


( 0.5, 


0.5, 


0.0) 


( 1.0, 0.5, 


0.5) 


( 1-0, 


0.5, -0.5) 


-0.00038 


8 


( 0.5, 


0.0, 


0.5) 


( 0.5, 1.0, 


0.5) 


( 1-0, 


0.5, 0.5) 


-0.00019 


9 


( 0.5, 


0.5, 


0.0) 


( 1.0, 1.0, 


0.0) 


( 1-0, 


0.5, 0.5) 


-0.00096 


10 


( 0.5, 


0.5, 


0.0) 


( 1.0, 1.0, 


0.0) 


( 1-5, 


1.5, 0.0) 


0.00120 


11 


( 1-0, 


0.0, 


0.0) 


( 1.0, 0.5, 


0.5) 


( 1-0, 


0.5, -0.5) 


0.00022 


12 


( 0.0, 


0.5, 


0.5) 


( 0.5, 0.0, 


0.5) 


(-0.5, 


-0.5, 1.0) 


-0.00043 


13 


( 1-0, 


0.0, 


0.0) 


( 1.0, 0.5, 


0.5) 


( 0.5, 


1.0, 0.5) 


-0.00005 


14 


( 0.5, 


-0.5, 


0.0) 


( 0.5, -0.5, 


-1.0) 


(-0.5, 


-1.0, -0.5) 


0.00005 


15 


( 0.5, 


0.5, 


0.0) 


( 1.0, 1.0, 


0.0) 


(1.0, 


0.0, 0.0) 


0.00032 
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FIG. 2: (Color online) Concentration dependence of the chemical interactions for the first five 
coordination shells. The lattice constant for each concentration is equal to the corresponding 
experimental value at room temperature. 
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FIG. 3: (Color online) Volume (lattice constant) dependence of the chemical interactions for the 
first five coordination shells at 17 at% Mn. The dashed line corresponds to the room-temperature 
value of the lattice constant for this composition. 
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FIG. 4: Left-bottom panel: Fourier transform of the SRO parameters obtained from Monte Carlo 
simulations for Cuo.83Mno.17 at a temperature 500K. Right-top panel: The atomic part of the 



intensity from a neutron-scattering experiment (Ref. 



. Dark regions correspond to maxima in 



the intensity; light regions, to minima. Wave vectors are given in reciprocal-lattice units (r. 1. u.). 
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FIG. 5: (Color online) Ordered structure DOqq. Light spheres represent the Mn atoms; dark, the 
Cu atoms. Smaller and larger spheres depict, respectively, positions at the corners and at the face 
centers of an elementary cubic cell. The elementary cell of the ordered structure is marked with a 
frame. 
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FIG. 6: Left-bottom panel: Fourier transform of the spin-spin correlation function obtained from 
Monte Carlo simulations for Cuo.83Mno.17 at a temperature 20K. Right-top panel: The magnetic 
part of the intensity from a neutron-scattering experiment (Ref 18|). Dark regions correspond to 
maxima in the intensity; light regions, to minima. Wave vectors are given in reciprocal lattice units 
(r. 1. u.). 
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FIG. 7: (Color online) Finite-size scaling of the SG susceptibility. 
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FIG. 8: (Color online) Finite-size scaling of the SG correlation length. Curves for different system 

sizes cross at Tc = 57 ± 5K. 
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FIG. 9: (Color online) Universal behavior of the SG correlation length. Best fit is achieved with 
V = 0.95. 
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FIG. 10: (Color online) Universal behavior of the SG susceptibility. Best fit is achieved with 
r/ = 0.25. 
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FIG. 11: (Color online) The extended scaling of the SG correlation length according to Eq. (fT7|) . 

The reduced temperature t is defined as t = T/Tc. 
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FIG. 12: (Color online) The extended scaling of the SG susceptibility according to Eq. (fT7|) . The 
reduced temperature t is defined as t = T /Tc. 
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